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Abstract 



t> ■ We show that for [/(l)-symmetric spacetimes on S* 3 x R a constant of mo- 

■ tion associated with the well known Geroch transformation, a functional 

OS 

■ K[hij, 7T*-'] , quadratic in gravitational momenta, is strictly positive in an open 

m ' 

subset of the set of all [/(l)-symmetric initial data, and therefore not weakly 

O ' 

zero. The Mixmaster initial data appear to be on the boundary of that set. We 
calculate the constant of motion perturbatively for the Mixmaster spacetime 

■ and find it to be proportional to the minisuperspace Hamiltonian to the first 
C$ • order in the Misner anisotropy variables, i.e. weakly zero. Assuming that K is 

exactly zero for the Mixmaster spacetime, we show that Geroch's transforma- 
tion, when applied to the Mixmaster spacetime, gives a new [/(l)-symmetric 
solution of the vacuum Einstein equations, globally defined on S 2 x S 1 x R, 
which is non-homogeneous and presumably exhibits Mixmaster-like compli- 
cated dynamical behavior. 
PACS number (s): 04.20.Cv, 04.20. Jb 
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I. INTRODUCTION 



The homogeneous cosmological models, apart from their possible cosmological relevance, 
provide us with a rich theoretical laboratory in which, it is hoped, many of the problems 
of General Relativity can be reduced to a manageable level of complexity One of the 
longstanding problems of Classical Relativity that benefited from such an approach is the 
problem of behavior of the solutions of Einstein's equations close to the "singularity", i.e., 
close to the boundary of the maximal globally hyperbolic development; a problem which 
is relevant to the issue of Strong Cosmic Censorship and possibly to the quantization of 
gravity program. Due to the complex nonlinear nature of the Einstein partial differential 
equations, an exact and full description of the asymptotic behavior seems to be unattainable, 
at least so in the forseeable future. Restricted to the spatially homogeneous cosmological 
models, however, the Einstein field equations, in the natural symmetry-adapted foliation, 
reduce to a system of ordinary differential equations which is much easier to analyze, and 
had been extensively analyzed in the past (see [|1||). Among the homogeneous models the 
most important special cases for our theoretical laboratory are the Kasner and Mixmaster 
models. The Kasner solution is explicitly known and is a very simple (generically) curvature 
singular solution. The Mixmaster model, despite several decades of intense effort, has not 
yet been explicitly solved (we do not have as detailed knowledge of its properties as we would 
want to), but an approximate description for the approach to the singularity has been found 
in terms of a discrete sequence of Kasner solutions f2|. Due to the stochastic properties 
of the associated discrete mapping it is hard to estimate the quality of the approximation 
analytically, but numerical simulations indicate that the discrete sequence approximates the 
exact solution the better the closer we come to the singularity j3]. 

The significance of the Mixmaster model for the problem of the asymptotic behavior of 
solutions of Einstein's equations stems largely from the long series of papers by Belinski, 
Khalatnikov and Lifshitz hereafter BKL. They tried to describe the approach to the 
singularity for non- symmetric spacetimes as a sequence of pointwise Mixmaster-like tran- 
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sitions between Kasner epochs. A satisfactory geometrical formulation of their method is 
still lacking, and it is still controversial whether their method can give any information 
about the global structure of the singularity ||. Nevertheless, it indicates that even lo- 
cally the dynamical behavior close to the singularity can be extremely complex; at least as 
complex as in the Mixmaster model (in fact as complex as in non-diagonal Bianchi IX and 
Bianchi VIII). Even though the asymptotic dynamics can be extremely complex for generic 
solutions, there exists a class of solutions whose asymptotic behavior is fairly simple: the 
asymptotic dynamics simplifies significantly for spacetimes for which, in a suitable foliation, 
the spatial derivative terms can be neglected asymptotically. Since the thus truncated Ein- 
stein's equations are exactly solvable we can extract all the asymptotic properties we need 
from the solution of the truncated equations, the Generalized Kasner Solution; a solution 
which evolves pointwise like the Kasner solution. The corresponding approximation method, 
usually called the Velocity-Dominated Approximation (VDA) |||7| or Strong Coupling Ex- 
pansion has been shown to be applicable to a large class of non-homogeneous spacetimes 
including Gowdy spacetimes Although the class of spacetimes for which the VDA is 
valid is infinite dimensional, it does not include some of the homogeneous models. Namely, 
the necessary conditions for the applicability of the VDA are not satisfied for the Bianchi 
VIII and IX spacetimes and, as BKL showed, are not satisfied for a generic spacetime as 
well. 

In order to test whether the Mixmaster solution (diagonal Bianchi IX), or the more 
general non-diagonal Bianchi VIII and IX, as claimed by BKL, are in some sense good 
asymptotic, pointwise approximations for the generic spacetimes, as the Kasner solution 
is for the velocity-dominated spacetimes, it is important to have some non-homogeneous 
solutions that are close to the Mixmaster in the gravitational phase space. If this is the case 
it is also important to have as detailed knowledge of the complex Mixmaster dynamics as 
possible. 

To that end in this paper we use the structure associated with the well known Geroch 
solution generating technique for spacetimes with one Killing field |§. Most importantly, 
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with Geroch's transformation there is associated a new non-trivial constant of motion, a 
spacetime observable which, we hope, could help us in the analysis of the asymptotic behavior 
of spacetimes. We treat the Mixmaster spacetime as a special case of [/(l)-symmetric 
spacetimes by using only one of the three Mixmaster Killing fields — one that generates a 
one-dimensional isometry group U(l). The action of that U(l) group on any Mixmaster 
homogeneous spacelike slice E, E ~ S* 3 , induces a Hopf fibration of E, i.e., makes it into 
the n — 1, U(l) principal fiber bundle over S 2 . We then calculate perturbatively the new 
constant of motion for the Mixmaster spacetime which, to the first order in the Misner 
anisotropy variables, turns out to be proportional to the minisuperspace Hamiltonian. This, 
barring an unlikely coincidence, indicates that the constant of motion is proportional to the 
minisuperspace Hamiltonian. If that is true, we can, using Geroch's transformation, obtain 
a new globally (in space) defined, non- homogeneous solution on S 2 x S 1 , with presumably 
Mixmaster-like complex asymptotic behavior. 

However, when calculated for the Taub-Nut special case, a spacetime that has 4 Killing 
fields, the constant of motion associated with the fourth Killing field, which does not exist 
in the general Mixmaster case, is strictly positive and therefore not proportional to the 
Hamiltonian constraint (which is zero for the solutions of the Einstein equations). It is 
worth emphasizing that there is no contradiction in the Mixmaster constant of motion being 
zero and the Taub-NUT constant being positive, since the above mentioned Taub-NUT 
constant of motion is not a restriction of the Mixmaster constant of motion to the Taub- 
NUT symmetry-class because a different Killing field is used for the Taub-NUT case. 

The outline of the paper is as follows: in section 2 we review the partial reduction of the 
Einstein equations in the ADM hamiltonian formulation for the [/(l)-symmetric spacetimes 
as well as the conditions for global applicability of Geroch's transformation. In section 
3 we describe the Mixmaster spacetime as a [/(l)-symmetric spacetime. In section 4 we 
apply Geroch's transformation to the Mixmaster and Taub-NUT spacetimes and find their 
respective constants of motion. In section 4 we give some concluding remarks. 
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II. [/(l)-SYMMETRIC SPACETIMES AND GEROCH'S TRANSFORMATION 



As shown by Geroch, the vacuum Einstein equations for the class of four dimensional 
spacetimes with at least one Killing field can be reduced, locally, to the three dimensional 
Einstein equations coupled to a harmonic map, with the two-dimensional hyperbolic space 
as the target manifold for the harmonic map. The isometry group, SL(2,R), of the target 
space becomes a symmetry of the reduced equations, i.e., maps solutions of the vacuum 
Einstein equations into new locally defined vacuum solutions with one Killing field. 

It was hoped that, with some control over the global structure of the spacetime, the 
transformation might be applicable globally. Indeed, for the globally hyperbolic spacetimes 
foliated by spacelike, compact, connected and orientable 3-surfaces (Cauchy surfaces) £, 
invariant with respect to the action of the isometry group, the conditions for generation 
of new globally (in space) defined vacuum solutions were found by Moncrief and Cameron 
|10| , |TT|| . They showed that the Einstein equations for U(l)-symmetric spacetimes on U(l) 
principal fiber bundles (circle bundles) Sxi? — > T,xR — where the base manifold £ ~ H/U(l) 
is a compact, connected and orientable two-dimensional manifold — can be reduced to a 2 + 1 
Einstein-Harmonic map system from £ x R to the Poincare half-plane (hyperbolic space) 
as the target space, provided one integrality condition is satisfied [|IC| , [n}| . For each base 
manifold £ there is a countable infinity of isomorphism classes of principal £/(l)-bundles 
B n over S, and the integrality condition depends only on the integer n characterizing the 
isomorphism class B n of the Cauchy surface E. For E = S 2 all classes were explicitly 
constructed by Quiros et. al. fl2| by imposing suitable identifications on the S 3 . The 
integer n characterizing the classes is equal to both the winding number and the Chern 
number of the bundle; the trivial bundle S 2 x S 1 ~ Bo, i.e., n = and S 3 ~ B\, i.e., n — 1, 
for example. 

The restriction to the U(l) isometry group, i.e., the restriction that the group orbits be 
closed, is not a restriction for the spacetimes with at least one dimensional isometry groups. 
The isometry group of a cosmological spacetime must be a compact Lie group — because 
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it is also the isometry group of the foliation surfaces E, which are compact Riemannian 
manifolds — and every compact Lie group has a U(l) subgroup. The only real restriction 
is that there exists a U(l) subgroup of the isometry group whose action on £ makes the 
foliation surfaces S into U(l) principal fiber bundles with compact and orientable base 
manifold S, i.e., that there are no fixed points for the action and no orbits twisting around 
each other. 

We shall now, mostly following ||10|| , give a brief review of Geroch's formalism. In order to 
express the conditions for the global applicability of Geroch's transformation as conditions on 
the gravitational initial data, we shall use the usual ADM hamiltonian formulation adapted 
to the class of spacetimes with one spacelike Killing field. Depending on convenience, we 
shall use index-free or typical-component (abstract-index) notation for tensor fields. The 
lower case greek letter (/x, v • • •) indices shall denote components in an arbitrary basis and 
will run from to 3. When dealing with coordinates (or tetrads) specially adapted to a 
family of spacelike surfaces we shall use the lower case latin indices, beginning with i, for 
the range from 1 to 3, and call such components the spacelike components. Likewise, when 
using further specialized coordinate systems (or tetrads), specially adapted to the group 
action, we shall use the lower case latin indices, beginning with a, for the range from 1 to 2. 

Let the U(l) isometry group be parametrized by the usual angle parameter a G [0, 2tt) 
and let the Killing field £ be just the derivative with respect to the group parameter a, i.e., 
£[/(p)] = ^/( e * Q p); where e ia p denotes the action of the U(l) group on the point p. By 
definition, the spacetime can be foliated by a family of spacelike hypersurfaces S t invariant 
with respect to the action of the isometry group, i.e., a family whose defining global time 
function t has zero Lie derivative with respect to the Killing vector field £ generating the 
isometry group action. The Killing field is then necessarily tangent to the foliation. 

As usual, using the unit normal vector field n M to S and the spacetime metric g^ u one 
can define a new tensor field, the three-metric 

h/a, = g^v + fiVA 9vo n x n a , (1) 



6 



whose restriction to the vectors tangent to £ generates a Riemannian metric on E. In a 
coordinate system generated by extending arbitrary coordinates on E using the flow of the 
normal vector n^, only spacelike components hij are non-zero, and both the three- metric 
h^ v and the foliation unit normal vector field n M are group invariant. 
The three-metric h^ u can be further decomposed as follows: 

h^u = A^ 1 ^ + A r]^, (2) 

using the norm A of the Killing field 

a = g^ec = K„ee, (3) 

and a one-form r)^ defined as 

V, = A" (4) 

This decomposition is ^-invariant and induces a new Riemannian two-dimensional metric 
g^v in the subspace of vectors tangent to £ and orthogonal to £. As a result one gets a 
^-invariant decomposition of the metric g^ v — > (n^, g^, rj^, A), and it is easy to see that 

= o, (5) 
e~9„ = 0, (6) 

ev, = i, (7) 
= o. (8) 

The last equation being a consequence of (|7|) and £^77 = 0. 

The bundle projection pullback *ir generates a one to one correspondence between, on 
the one side, the ^-invariant covariant tensor fields onS ~ B n all of whose contractions with 
£ vanish, and on the other, the covariant tensor fields on the base £ [[OJ, We can, therefore, 
treat the fields A and g a b as tensor fields on £ (and Ti^ clS cl field onExfl). The one-form rj, 
however, does not have vanishing contraction with £ and cannot be treated as a field on E, 
which imposes the only global restriction to the reduction of dynamics from E to E. Even 



though f] does not have vanishing contraction with £, its exterior derivative (the only form in 



which rj will enter the equations of motion) does, and as shown in |TJ[ , for a closed two-form 
$ on E there will exist a one-form 77 on E ~ B n , satisfying ^rj^ = 1 and drj = *7r($), if and 
only if 



f $ = 2vrn. (9) 



To effectively use the symmetry of the problem in concrete calculations we shall be 
doing in this paper, we need an isometry-adapted atlas on the spacetime manifold. Using 
an atlas on the base manifold E, with coordinates in a representative chart labeled by {x a }, 
a = 1,2 one can construct, with the help of the bundle structure on E, a ^-invariant atlas 
on the initial foliation surface E. Let the coordinates in a representative chart be labeled 
{x 1 } = {x a ,x 3 }, i = 1..3, with the third coordinate normalized so that d/dx 3 = £. Now, 
one can choose a ^-invariant time-evolution vector field t M and lift the coordinates from the 
initial leaf E, along the flow generated by £ M , to all other foliation surfaces, which are labeled 
by the flow parameter t used as the zeroth coordinate. Thus, one has created an atlas on 
the spacetime manifold in whose every chart, = {t,x a ,x 3 }, // = 0..3, the action of 

the U(l) group is given by e ia (t,x a ,x 3 ) = (t,x a ,x 3 + a), and the bundle projection % by 
7i : (t, x a , x 3 ) h-> (t, x a ). 

In any such chart we can define a one-form 

(3 = 7]- dx 3 , (10) 

which carries all important dynamical information contained in r] and whose exterior deriva- 
tive, d/3 = drj, and time derivative, $ = C t (3 = r), are globally defined and chart independent 
differential forms. The two-form d(3 can be identified with the (exact, globally-defined) 
two- form $ on the base manifold E and is, in an x a chart, expressible as 

d/3 = drj = r dx 1 A dx 2 , (11) 

where r is a scalar density whose value in the same chart, 



r = e ab (3 a , b , 



(12) 



can be calculated using the contravariant antisymmetric tensor (density) e ab = —e ba , e 12 
associated with the chart. 

In the new ?7(l)-adapted variables the isometry of the spacetime metric, 

d s 2 = g„ v dx^dx u = -N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) 

= X- 1 \-Ndt 2 + g ab (dx a + N a dt)(dx b + N b dt) 



+A dx 6 + [3 a dx a + N a (3 a dt 



(13) 



is equivalent to all the component functions being independent of x 3 , as the coordinate 
one-forms {dt, dx 1 } are already £- invariant. 

Using the U(l) isometry, the ADM action / on the E x R, 



I = J dtj^dx 3 {ir ij h ijyt — NH — WHi} , 



(14) 



with 



U = h- 1 ' 2 [vr%- - l/2(7ri) 2 ] - h l ' 2 R{h), 
Hi = -2/i 1 /2 V .( /l -i/2 7r i )) 



(15) 
(16) 



reduces to an equivalent action / on £ x R, 



where 



I = 2tt J dt jf dx 2 {n ab g ab , t + p\ t + e a /3 a , t - NH - N a H a + f3 e a a ) , (17) 



H = \-V 2 H = g- 1/2 



2 i 1 \2 2 , 1 ~ ah 



n ab n ab - (n a a ) + g A V + ^^e a e 



+ ~9 1/2 



-R(g) + ^9 ab \a\ b + ^~9 ac ~9 bd (d(3) ab (d(3) dc 



H a = H a = -2~g l ' 2 V b (g- 1/2 t b ) + p\, a + e b (d(3) ab , 
N = A 1/2 iV, 
A> = N a , 



(18) 
(19) 
(20) 
(21) 



A> = ViN\ (22) 

e a = 2\r ll n ia , (24) 
p = 2^'^^- - A ; ;i ;, 7/, r (25) 

Variation of the action J with respect to flo gives a constraint: 

e a a = ^ div e = ^ d(*e) = 0; (26) 

where we have used e = g -1 / 2 e to denote vector field associated to the vector density e, 
and * to denote the Hodge star operator of the Riemannian metric gon E. This constraint 
can be easily solved since, according to the Hodge theorem, every closed form on E, and in 
our case this is the form *e, can be uniquely decomposed to a sum of one exact and one 
harmonic form. Explicitly written in an arbitrary chart on E, the decomposition gives 

e b e ba = u, a + h a , (27) 

with uj the scalar field (defined up to a constant) and h a the components of the unique 
harmonic form, both defined globally on E. By a slight abuse of notation, we have used 
e ab = —e ba , £i2 — 1, to denote the covariant antisymmetric tensor (density) associated with 
the chart. From the equations of motion obtained by varying the action with respect to f3 a , 
we can obtain 

c^,a + K = \N~g- 1/2 \ 2 r + N\iu, b + h)] , (28) 

which, because of the uniqueness of Hodge decomposition and the fact that the one form 
on the right-hand side is exact, forces the harmonic form h to be zero, i.e., h to be time 
independent. Hereafter in this paper we shall be interested only in the special case of E ~ S 2 
which has no non-zero harmonic one-forms, and shall therefore put h — 0. 

With the decomposition of e a included into the action I, the one-form f3 a appears only in 
the form e ah fi a ^ = r and the action reduces to an equivalent 2+1 Einstein-Harmonic action 
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ruj.t - 



NH - N a H a } , 



(29) 



where 



H = g 1 7T 7Y ab - [7TJ + -G Fa^B 



+ g 1 ' 2 \-R(g) + -G AB X A a X B b g ab , 
H a = -2~g l ' 2 S/ b (g- 1/2 K b a ) + p\a + ruj A 
X A = {X,oo} } P A = {p,r}, A = 1,2, 



(31) 



(30) 



(32) 



and 



A" 2 d\ 2 + dw 2 



(33) 



is the metric of constant negative curvature on the Poincare half-plane {A > 0,u}, the 
target space for the harmonic variables A and u. This form of the action makes evident the 
SL(2,R) symmetry of the equations of motion, SL(2,R) being the isometry group of the 
Poincare half-plane, and enables us to apply the Geroch transformation to the harmonic 
variables A and u to generate new solutions of the vacuum 3+1 Einstein equations from the 
known ones. 

Starting with a [/(l)-symmetric four-dimensional solution of the vacuum Einstein equa- 
tions on Rx B n expressed in the form ([13|), we can calculate the canonical variables (A, r, p) 
locally from the metric components and their velocities. A is the norm of the Killing vector 
field, r is given by fll2"D and 



Only the canonical variable u is not a local function of the metric components and their 
velocities (or momenta), but it can be calculated by performing a line integral in the three- 
dimensional space i?xE: 




(34) 




(35) 
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In order to do so, we need the spatial and temporal derivatives of uj. Without using the 
equations of motion in the Einstein-Harmonic form the spatial derivatives can be calculated 
from (|24| ) and ( p7[ ) and the time derivative from fl28|) giving: 



uj >t = N~g~ l/2 \ 2 r - N-'N^W, ~g hc [5 c + c(t), 



(36) 
(37) 



where c(t) is an arbitrary function of time. This is to be expected, since uj is defined only up 
to a constant on each E, i.e., given anwa whole class of functions uj + f(t) is associated to 
each four dimensional spacetime. So, it would appear that there is a gauge indeterminacy 
for the evolution of uj. Fortunately, the new evolution equations completely eliminate that 
gauge freedom for uj and, given initial data for X,uj,p and r, determine a unique solution. 
Varying the action J with respect to r, we find that the class representative for uj singled 
out by the new equations of motion is the one obtained by putting c(t) = in (|3"7|). The 
class representative is not completely determined — a constant (in time and space) can be 
added to u — which corresponds to freedom of choice of the reference point xq in the initial 



E from which the spatial line integration in ( 3q) starts. 
The action of an element, 



y c d j 



ab — cd 



(38) 



of SL(2, R) on the canonical variables, given by: 

A 



A' 



UJ — > UJ 



p — > p 



c 2 (uj 2 + \ 2 ) + 2cduj + d 2 ' 
ac(uj 2 + A 2 ) + (ad + bc)uj + bd 
c 2 (uj 2 + A 2 ) + Icduj + d 2 ' 
p [c 2 (uj 2 - A 2 ) + Icduj + d 2 ] - r [2\{cd + c 2 uj)} 

c 2 {uj 2 + A 2 ) + 2cduj + d 2 
p(2c 2 Xuj + 2cd\) + r \d 2 + c 2 {uj 2 - A 2 ) + 2cduj 



(39) 

(40) 

(41) 
(42) 



then transforms our old solution of the 2+1 Einstein-Harmonic equations on R x E into 
a new solution of the same equations. To lift that solution to a new solution of the 3+1 
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vacuum Einstein equations on R x B n /, all we need to do is use new A' and integrate ((TTJ) to 
find the new one-form vj '. As mentioned previously, this will be possible if and only if the 
integrality condition (|9]), which now becomes 

r' = 2im', (43) 



s 

is satisfied at all times. 

The above integral is a constant of motion of the Einstein-Harmonic equations and the 
integrality condition has to be enforced only at one time. Even better, the symmetry-action 
of SL(2, R) results in three conserved quantities, each associated with one of the three Killing 
fields of the Poincare half-plane. These three constants of motion are: 

A = J_(2ujr + p), (44) 
# = /.r, (45) 
C = J^[r(\ 2 -u 2 ) -pu], (46) 

and represent a momentum map of the SL(2, R) action. 

B is not interesting as a constant of motion since its value does not depend on the initial 
data; it is just the integrality condition dictated by the topology of S. A and C are not 
proper constants of motion associated to the spacetime, either, because their values depend 
on the constant added to ou, which is not observable in the spacetime. Nevertheless, the 



SL(2, _R)-invariant constant of motion found by Geroch |13 | 



K\h l3 y ] \ = A 2 + ABC. (47) 

has the desired properties. It is easy to see that the value of K, not only does not depend 
on the unobservable constant added to u, but remains the same even if we add an arbitrary 
function of time to ou, which means that it is a proper constant of motion of the spacetime, 
i.e., an observable. This also means that when calculating u to be used for evaluating K, 
the first term in fl3"5D, which is only a function of time, can be dropped. That guarantees 
that the constant of motion K depends only on the initial data on a single Cauchy surface, 
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i.e., that it is a functional of the standard gravitational variables (h^, 7T lJ ) that is local in 
time. 

In order to calculate K[hij, 7r y ], given the ADM initial data (hij, 7r y ), one has to express 
A, p, r and uj ja in terms of the components of the ADM data in a 17(1)- adapted coordinate 
chart: 

A = h 33 , (48) 
p = 2\-\^h i3 h j3 - X-Whij, (49) 
r=(A- 1 / i 3i) 2 -(A- 1 / i32 ) ii , (50) 
oj >a = 2e ba n bi h i3 , (51) 

and evaluate u by integrating uo A along a curve in the two-sphere. The constants of motion 
A, B and C are then calculated by evaluating the surface integrals. Loosely speaking, K is 
a quadratic functional of the ADM momenta 7r lJ . 

An important property of K is that it controls the global applicability of the Geroch 
transformation. As shown by Moncrief [ITS], any solution on a nontrivial bundle B n , n ^ 
0, can be transformed to a globally defined solution on any other nontrivial bundle B n i, 
n' 7^ 0; no restrictions exist. For a more interesting case of transforming a solution from a 
nontrivial bundle to a global solution on the trivial bundle there is a restriction, however. 
The transformation can generate a globally defined solution if and only if K > 0. 

In the next section we shall establish some of the properties of K and try to calculate it 
for the Mixmaster spacetime. 

III. MIXMASTER MODEL AS A SPECIAL CASE OF THE [/(l)-SYMMETRIC 

SPACETIME 

In a spatially homogeneous spacetime, by definition, the orbits of the isometry group 
(or some subgroup if there are timelike Killing vector fields) are three-dimensional space- 
like hypersurfaces that foliate the spacetime. If the isometry group has a subgroup that 
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acts simply transitively on the orbits, the spacetime can be classified into one of the nine 
Bianchi classes, according to the Bianchi class of the Lie algebra of the corresponding simply 
transitive subgroup. If the isometry group has no simply transitive subgroup, the isometry 
group must be at least four- dimensional and the spacetime must be locally isometric to 
the interior of the Schwarzschild solution. These solutions, found by Kantowski and Sachs 
16fl , if restricted to compact spatial sections must have the spatial sections diffeomorphic to 



5 2 x S 1 and the four dimensional group with transitive action must be SU(2) x U(l). 
The Mixmaster cosmological model, by definition, belongs to the Bianchi IX class, has 

5 3 spatial topology and SU(2) as the isometry group, and the spatial metric is diagonal in 
the left- or right-invariant basis of the isometry group. 

Let us begin with the brief analysis of the structure of the spatial sections. The three- 
sphere, the group manifold of SU(2), can be represented as a hypersurface in C 2 ~ R 4 
S 3 = {(^1,^2), \zi\ 2 + I ^2 1 2 = 1}- Writing z\ = r\e ai and z 2 = ne" 2 , with > and a« 
from any interval of width 2ir, we see that the only condition the points lying on S 3 must 
satisfy is (rx) 2 + (r 2 ) 2 = 1, which guarantees that there exist unique 9 G [0,7r] such that 
r 1 = cos6'/2 and r 2 = sin6 l /2. Thus, one can introduce {9,ai,a 2 } as coordinates on S 3 , 
and can think of S 3 as a singular one-parameter family of two-tori (a\,a 2 ), 9 being the 
parameter of the family. The tori for 9 = 0, 7r collapse to S 1 , making the coordinate system 
singular there. 

Recall that we need a coordinate system in which one coordinate vector field would 
have closed orbits, i.e., generate a U(l) group action, and that group action would have to 
be such that it gives S 3 a principal fiber bundle structure over some two-dimensional E. 
The most useful coordinates on S 3 , for that purpose, are the standard Euler angle coordi- 
nates {9,<j),if}}, which can be obtained from {9, ai,a 2 } coordinates by a simple coordinate 
transformation: 



(j) = OL\ + «2, 1p = d\ — OL2- 

The Euler coordinates, just like the {9, ai,a 2 } have a singularity at 9 
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(52) 



= 0, 7r and when we 



let {9, 0, ip} e {(0, tt), [0, 2vr), [0, 4vr)} or {0, 0, ^} G {(0, vr), [0, 4tt), [0, 2tt)} they cover all of 
the three-sphere (except the two singular circles) exactly once. The period of both and 
ip coordinates is 4tt, and if we let both and ip extend all the way to 4-7T, the three-sphere 
would be covered twice. 

The coordinate vector fields 8$ and 8^ are globally defined on S 3 despite the coordinate 
singularities in 9 = 0, tc and they generate two U(l) group actions, which each yield a Hopf 
fibration of S 3 , i.e., they make S 3 into a principal fiber bundle over S 2 with bundle projection 
for the ^-bundle, for example, given by: 

n^:S 3 ^S 2 , (#,0,^)^(M)- (53) 

Since these two bundle structures are equivalent, we shall use ^-bundle for the Mixmaster 
spacetime and, taking -0/2 = x 3 , define 

£ = 2 fy, (54) 

to be the Killing field we shall use to apply Geroch's transformation to the Mixmaster so- 
lution. The coordinate system {x x ,x 2 ,x 3 } = {9, 0, ip/2} is appropriately adapted to the 
bundle structure, as described in the previous section. In a spacetime with several Killing 
fields it is important to remember which Killing field was employed for Geroch's transfor- 
mation, so we shall denote the constant of motion ([17]) associated with the Killling vector 
2 8^ by K^. 

The set of standard, globally-defined (and analytic) left-invariant vector fields X iy which 
generate a right-action on 577(2) |17|| , written here together with their dual one-forms lj 1 : 

X\ = cos ip do + esc 9 sin ip 8$ — cot 9 sin i[> 8^ 

X 2 = — sin ip dg + esc 9 cos ip 8$ — cot 9 cos ip 8^ 

X-i = 8,1, 

a) 1 = cos ip d9 + sin 9 sin ip d<fi 

Cj 2 = — sin ip d9 + sin 9 cos ip d(f) 

Cj 3 = dip + cos 9 d(p 
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ci x k 



(56) 



Xi, Xj 

dCj k = -1/2 C% u* Au j 

C\j = Cijk, £123 = 1 

contains d^, so, we shall chose them as the complete set of Killing fields for the Mixmaster 
spacetime. We shall also need a right-invariant one-form basis a 1 to explicitly express the 
isometry of the metric. For completeness we give them with their dual vector fields % and 
some of their properties: 



Y 1 



a 1 



cos (ft de + esc 9 sin (ft — cot 9 sin (ft 
sin (ft de — esc 9 cos (ft + cot 6* cos (ft 

cos (ft d9 + sin 6* sin dip 
sin (ft d9 + sin 6* cos oh/> 
d0 + cos 9 dip 



(57) 



y y 



da k = 1/2 C% a i A a j 



Xi, Yj 



(58) 



: 0, C x°j = 

The Mixmaster metric, invariant with respect to the right action of SU(2), can now be 
written as: 



ds * = ,^2 dt 2 + A 2^2 + S 2^2 )2 + £2^3^ 



(59) 



where AT, A, _B and C are functions of time only, and when expressed in the usual Misner 
anisotropy variables Q, j3 + and 



5 



(60) 
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IV. THE CONSTANTS OF MOTION AND TRANSFORMED SPACETIMES 



As the first step in calculating the constants of motion for the Mixmaster spacetime, we 
need to calculate the harmonic variables. The function A and one-form (3 can be found simply 
by expanding the metric (0). The functions r and p can be straightforwardly calculated 
using fll^D and (|51D, which requires only differentiation. The function lu, on the other hand, 
is hard to calculate since that involves a line integration in the (#,</>) surface with 

A = A 2 sin 2 sin 2 6 + B 2 cos 2 sin 2 9 + C 2 cos 2 6 (61) 

in the denominator. We have not been able to evaluate u in closed form, which prevented us 
from calculating the constant of motion exactly. But, even if one succeeded in putting 
uj into a complicated closed form it would probably be of little practical importance, since 
we would still have to evaluate the surface integrals, which seem extremely difficult. 

However, the integrals are simple enough when A = B = C, i.e., when the anisotropy 
variables j3 + = (3- = 0, which allowed us to calculate perturbatively around the isotropic 
solution, to the first order in the anisotropy variables. This calculation is tedious and rather 
unilluminating, so we shall state only the main result: 

To the first order in the anisotropy variables (/3 + , f3J), the constant of motion for the 
Mixmaster spacetime is proportional to the minisuperspace Hamiltonian, 

H = 1/12 Ne 3n [-pi +p 2 + + p 2 _ + e~ 4n V((3 + , /?_)] , (62) 
V(J3+, /?_) = 3 {e- 8(3+ - 4e" 2/3+ cosh(2 v / 3/3-) + 2e 4/3+ [cosh(4^3/3_) - l] } . (63) 

In particular, for all Mixmaster solutions whose trajectories in the anisotropy plane at least 
once pass through the origin, j3 + = (3_ = 0, the constant of motion = exactly. 

Is K^p exactly proportional to the minisuperspace Hamiltonian? We do not know, but 
it seems to us that, considering the complexity of the calculations involved, it would be an 
extraordinary coincidence if the first two terms in the expansion of K coincided with those 
of H, without K and H being essentially equal (proportional). In addition, if all Mixmaster 
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trajectories passed through the origin the value of would have to be exactly equal to zero 
for all Mixmaster initial data. Regrettably, it is not known how large is the set of Mixmaster 
trajectories passing through the origin. 

Does this indicate that the constant of motion K is proportional to the hamiltonian 
constraint, i.e., uninteresting and trivial, for all [/(l)-symmetric spacetimes on S" 3 ? Is it 
weakly zero? Not at all! To see that, we shall calculate K for the Taub-NUT spacetime, 
using the additional Killing field the general Mixmaster spacetime does not have. 

The Taub-NUT spacetime is a special case of the Mixmaster spacetime, obtained by 



putting A = B in (p9|) . Its equations of motion can be explicitly solved, and give the 



solution (written here in the right-invariant form, slightly different form [18|, where it is 
given in the left-invariant form): 

ds 2 = -U- l dt 2 + {f + I 2 ) [(a 1 ) 2 + (a 2 ) 2 ] + (2/)2f/(a 3 ) 2 , 

_ -t 2 + 2mt + l 2 ( 64 ) 

Due to the additional restriction (A = B), the Taub-NUT spacetime has four Killing 
fields; the three left-invariant vector fields inherited from the Mixmaster plus one of the 
right-invariant vector fields, namely the 8$. Note that 8$, being a right-invariant vector 
field, commutes with other three Killing fields. This additional Killing field, as mentioned 
in the previous section, also generates a U(l) group action which induces a Hopf fibration 
of S 3 . We can use it (actually 2 8$ because of the normalization convention we use) to 
calculate a new constant of motion, which we shall call K$. Owing to the simplicity of the 
Taub-NUT solution all angular integrals can now be easily calculated and the constant of 
motion turns out to be: 

^[Taub] = 16(8tt) 2 (2/) 2 (/ 2 + m 2 ). (65) 

As we can see, -fQ,[Taub] is a strictly positive, non-trivial function of the initial data (I 
and m) and as such cannot be proportional to the Hamiltonian constraint, which is zero for 
all solutions of the Einstein equations. More generally, it guarantees that is not weakly 
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zero. This implies that there exists a neighbourhood O of the Taub-NUT initial data, open 
in the set of all ?7(l)-symmetric initial data on S 3 with as a Killing field, in which the 
constant of motion is strictly positive. That constant of motion is local in time as a 
functional of the usual gravitational initial data — the three-metric and its momentum — 
and can be written down essentially explicitly. It seems to be the only explicitly known 
observable for such a large family of Einstein spacetimes. 

For the left-invariant Mixmaster metric — which, as opposed to its isometric right- 
invariant counterpart (|59|) , has 8$ but not as Killing field — the K$ constant of motion can 
be calculated. (The left-invariant Mixmaster metric can be obtained from (|59|) by making 
A, B and C appropriate functions of the angles or exchanging as for uu-s.) Due to isometry 
of the left- and right-invariant spacetimes, it must be that [Mixmaster] = [Mixmaster] . 
Assuming that [Mixmaster] is indeed proportional to the Hamiltonian constraint, it fol- 
lows that the Mixmaster spacetime is the special case for which the [/( regenerated constant 
of motion is (weakly) identically zero. 

The positivity of fQ,[Taub] also implies that Geroch's transformation can yield a new 
globally defined spacetime on S 2 x S 1 when applied to Taub-NUT. Since the harmonic 
variables can be easily calculated in this case, the transformation can be carried out explicitly 
and we find that the new solution on the trivial bundle turns out to be the Kantowski-Sachs 
solution. This had already been established by Geroch himself ||13|| , but only locally. 

Now a word about the possibility of generating new globally defined solutions of the vac- 
uum Einstein equations by applying Geroch's transformation to the Mixmaster spacetime. 



Moncrief's analysis of the global applicability of Geroch's transformation |T5] guarantees 
that it is possible to transform a solution from S 3 , as the bundle with winding number 
n — 1, to a new solution on the trivial bundle S 2 x S 1 if K > 0. We showed that for all the 
Mixmaster solutions that pass through the f3 + = /3_ = point in the anisotropy plane K 
is exactly zero; the condition for global applicability of Geroch's transformation is therefore 
satisfied and we can obtain a new solution on the trivial bundle. Likewise, if K is exactly 
zero for all Mixmaster spacetimes, as suggested by our perturbative result, a new solution 
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globally defined on the trivial bundle can be obtained from any Mixmaster solution. 

What properties does that new solution have? It is not a homogeneous solution, for the 
following reason. If it were homogeneous it would have to be the Kantowski-Sachs solution, 
which is the only homogeneous solution on S 2 x S 1 . But the Kantowski-Sachs constant of 
motion K associated with the Killing vector field tangent to the S 1 factor is strictly positive. 
It is positive because it must be equal to the Taub-NUT constant K, since K is invariant 
under the action of Geroch's transformation and the Kantowski-Sachs solution is just the 
transformed Taub-NUT solution. The same invariance requires that the Kantowski-Sachs K, 
a strictly positive quantity, be equal to the Mixmaster constant K, which is identically zero, 
leading to a contradiction. So, the new solution on the trivial bundle is not homogeneous. 

The asymptotic behavior of the new solution requires additional analysis; for now let 
us just say that it is probably oscillatory in the BKL sense. The BKL condition for the 
existence of the oscillatory behavior is an "open" condition, i.e., some quantity has to be 
different from zero. It is satisfied by the Mixmaster solution, hence it has to be satisfied 
by the transformed solution obtained by using an element of SL(2, R) sufficiently close to 
unity. Whether the element of SL(2,R) that gives the new globally defined solution on the 
trivial bundle violates the BKL condition or not, will be the subject of future study. 

V. CONCLUSIONS 

The constant of motion K we calculated in the previous section is a second degree poly- 
nomial in the minisuperspace momenta and, to the first order in the Misner 
anisotropy variables it is proportional to the minisuperspace Hamiltonian. We 
believe that it is reasonable to conjecture that K is exactly proportional to the minisu- 
perspace Hamiltonian. The coincidental agreement of the first two terms in the expansion 
of K and H, considering the complexity of the expression for the constant of motion K, 
is highly unlikely unless K and H are exactly proportional. It seems that the Mixmaster 
dynamical system has once more defied the attempt to find a non-trivial constant of motion 
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that would help us to gain more detailed knowledge of its behavior. Of course, this might 
just be a consequence of the simple, but not yet established, fact that there are no non- 
trivial constants of motion and that a generic Mixmaster trajectory fills densely an open 
subset of the constraint surface in the minisuperspace. Indeed, Kuchaf |B| showed that 
for non-homogeneous spacetimes there are no (weakly) non-zero constants of motion that 



are linear functionals of the gravitational momenta. Presently K. Schleich and D. Witt j2G 
are trying to prove that any constant of motion for the Mixmaster dynamical system that 
is quadratic in minisuperspace momenta — which includes our K — has to be proportional 
to the minisuperspace Hamiltonian, i.e., has to be trivial. Our perturbative result for K 
supports that. 

More importantly, the probable weak vanishing of K for the Mixmaster spacetime does 
not mean that K is weakly zero for all [/(l)-symrnetric spacetimes. By calculating explicitly 
the constant of motion for the Taub-NUT spacetime, we showed that around the Taub- 
NUT intial data, there exists an open set in the set of all [/(l)-symmetric initial data for 
which the constant of motion K is a non-trivial spatially non-local, but temporally local, 
strictly positive functional of ADM canonical data. That constant of motion could be used 
as an observable in a t/(l)-symmetric toy model for the future quantum theory of gravity. 
Classically, the value of this constant of motion for a particular [/(l)-symmetric spacetime 
might signal whether the spacetime is asymptotically velocity-dominated or not. The value 
of K for the Taub-NUT spacetime is positive and the spacetime is velocity-dominated, 
on the other hand the value of K for the Mixmaster spacetime, which is not a velocity- 
dominated spacetime, is presumably zero, which might signal that the K = surface in the 
set of £7(l)-symmetric initial data is the boundary between the simple velocity dominated 
behavior and the much more complicated oscillatory BKL-like behavior. The Mixmaster 
spacetime would then be exactly at the boundary of the set of velocity-dominated spacetimes. 
Likewise, an interesting and closely related issue we would like to resolve is: Does Geroch's 
transformation, which does not change the value of K, "conserve" the velocity-dominated 
behavior, i.e., are new solutions obtained from velocity-dominated solutions also velocity- 
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dominated? 
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